AC induced damping of a fluxon in long Josephson junction 
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We present a theoretical and experimental study of Josephson vortex (fluxon) moving in the presence 
of spatially homogeneous dc and ac bias currents. By mapping this problem to the problem of 
calculating the current-voltage characteristic of a small Josephson junction, we derive the dependence 
of the average fluxon velocity u on the dc bias current 7. In particular we find that the low frequency 
ac bias current results in an additional nonlinear damping of fluxon motion. Such ac induced 
damping crucially depends on the intrinsic damping parameter a and increases drastically as a is 
reduced. We find a good agreement of the analysis with both the direct numerical simulations and 
the experimentally measured 7(m) curves of a long annular Josephson junction with one trapped 
fluxon. 
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I. INTRODUCTION 

Solitary waves have attracted a great attention for 
many years because they cajjij-|be excited and observed 
in diverse nonhnear systems. El'H Well-known examples of 
such solitary waves are optical solitons, domain walls, 
kinks and magnetic fluxons, just to name a few. 

The latter case, a fluxon carrying the magnetic flux 
quantum $0, has a form of a Josephson vortex mov- 
ing along a long Josephson junction (LJJ). This system 
presents an additional interest because the fluxon motion 
can be easily contrcjl|led by externally applied magnetic 
field and dc bias.cTEl It is also well known that under 
these conditions the fluxon motion can be mapped to the 
classical mechanics of a macroscopic (relativistic^ particle 
and such correspondence is practically exact .Bia 

It was shown theoretically and experimentally that in 
the presence of both dc and ac bias the fluxon motion 
becomes mO|Cje,cmnplex and it displays various fascinating 
phenomena.BljilJ Thus, as the frequency lu of ac bias is 



where ujp is the plasrua-firequency of 



large, i.e. lu 

a LJJ, the fluxon radiates plasma waves.crOlZl When the 
frequency u is small (w <^ ujp) but inhomogeneous (static 
or dynamic) potential is present, the fluxon motion can 
be locked to the external ac bias. In this case the vertical 
steps ( "Shapiro steps" ) appear in the dependence of the 
average fluxon vdodty u on the dc bias 7 (analog of 
the /- V curve). ffliJ Moreover, in the same frequency 
region the_resonant escape of a pinned fluxon has been 
observed. E2l 

However, the fluxon motion in the presence of both 
homogeneously applied dc and ac bias for the frequency 
uj <^ ujp, has not been investigated neither theoretically 
nor experimentally. In this Report we present a theoret- 
ical and experimental study of this limit. We show that, 
unexpectedly, the influence of a small ac bias is rather 
large, and an additional ac induced damping increases as 



the standard damping characterized by parameter a re- 
duces. 



II. EQUATION OF MOTION AND 
THEORETICAL ANALYSIS 

The dynamics of a fluxon in a long Josephson junc- 
tion is described by a perturbed, sin-Gordon equation for 
the Josephson phase (^(x, i).0O For a particular case of 
homogeneous dc (7) and ac (e sin cut) bias currents this 
equation has the forr 



(ftt — sin Lp = a(ft — 7 — £ sin ujt. 



(1) 



Here, the units of time and coordinate are correspond- 
ingly, the inverse plasma frequency 1 /ujp and the Joseph- 
son penetration length Aj. The dimensionless damping 
parameter a is determined by the direct losses that are 
due to the presence of quasiparticle (dissipative) current. 

In the absence of perturbation (rhs of Eq. (j^) ) the 
solution of Eq. (^) is a a fluxon moving with velocity u. In 
a general case we use the standard transformation of the 
Eq. (|^) to the fluxon reference frame, where fluxon moves 
with an average velocity u. New variables are defined as 

X -ut t-ux r 

p p 

The Eq. (|^) rewritten in the new variables is 

a an ( r + 

^ii - Vrr ~ sinyj = -^"^^ + - 7 - esm — - — 

(2) 

Note here, that the fiuxon, in its reference frame, "feels" 
the periodic coordinate and time dependent potential 
given by the last term of Eq. (0) . 
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We are interested in the regime of small frequencies of 
ac bias current, namely w <C 1 in the dimensionless units. 
In this limit we can neglect rthe interaction of the fluxon 
with the plasma oscillationsflu The potential energy iJpot 
of the fluxon is: 



+00 



- 7'/' + (1 - cos if) - eip sm Ilo — - — 1 



(3) 

Introducing the fluxon solution in the form r) — 
(po{^ — X{t)), where (Pq{^) is a well known Josephson 
vortex solution we derive the equation of motion for the 
coordinate X{t) of the fluxon (in the reference frame): 



•• a ■ au IT 

X + -X H = -7 - 
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2/3 
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(4) 



Thus, the main result can be seen from this equation: 
the applied ac bias current (last term in Eq. (||) ) excites 
fluxon librations of the frequency ^ and, in turn, the 
interaction of such librations with the ac bias current 
leads to the additional nonlinear damping of a fluxon 
motion. The interesting feature of the Eq. (Q) is that 
it can be mapped to the problem of calculating the I-V 
characteristic (IVC) of the small Josephson junction with 
a specific damping, in general, different from a. To see a 
precise mapping we introduce new variables: 



y 



-X. 



To 



UJ 



and obtain the equation 
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where ic{u) is a "critical current" of the small junction 
which depends on the velocity u of the fluxon from our 
problem and is given by: 



ic{u) 



■Kuj3e 



Auj cosh 
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(6) 



The 7-M characteristics of moving fluxon can be found 
from this equation in the following form: 



Aau 
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(7) 



where i{v) is the nonlinear current- voltage characteristics 
of a small Josephson junction with the critical current 
ic{u) and the damping a/uj. In Eq. (|^), the first term 



is the usual 7-u characteristic for a fluxon moving along 
LJJ and the second term corresponds to the correction 
due to the presence of ac bias current. 

The analytical expression for i{v) dependence for small 
JJ at arbitrary values of damping and critical cuxrent 
is not known. Using numerically simulated i(z;)E3 and 
Eq. (|^), we obtain the "f[u) dependence for various values 
c'Cof e. Fig. shows 7~M characteristics calculated in this 
way for e between zero and 0.7. 
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FIG. 1. The 7-u curves obtained by using Eq. (|7|) for vari- 
ous values of the ac bias amplitudes: e = 0, 0.1,0.3, 0.5, 0.7. 
The damping parameter a and ac frequency ui are, corre- 
spondingly, 0.04 and 0.03. The inset shows the normalized dc 
bias change A to the applied ac bias e = 0.1 at uq =0.3, as 
a function of the damping parameter a. 

Let us analyze some limits. For small fluxon velocity, 
namely u <C max{a, w}/e, 7(11) looks like a straight line, 
but its slope is steeper than in the autonomous case e = 
(see Fig. I): 



7 ^ 
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1 + - — 



f3' 



(8) 



In the opposite case u ^ max(a,ct;)/e, the 7- u depen- 
dence just shifts by the constant value e < 1 from the 
autonomous curve. Thus, ac bias induced damping is 
highly nonlinear, and the 7- u dependence displays an 
inflection point for moderate values of u (see Fig. 1). 

In the region of small frequencies w ~ a as the damping 
coefficient a decreases, the response of the fluxon to the 
ac bias current (which is proportional to 1/a^) drastically 
increases [see Eq. (^)]. In order to analyze this effect in 
more precise terms we calculate the change of dc bias 
^ = as a function of the damping parameter 

a, for particular values of fluxon velocity uq = 0.3 and ac 
bias e = 0.1, where 70(1*0) is the value of dc bias current 
in the absence of ac bias. The results are presented in 
the inset of Fig. |l]. 
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III. NUMERICAL SIMULATIONS 

In order to verify the analysis presented above and pre- 
dictions made we performed direct numerical integration 
of the PDE (|l|). The periodic boundary conditions 

^{e,t)^ipiO,t) + 2TT; (9) 
^o.{i,t) ^ ip,iO,t), (10) 

were imposed and the normalized length of the LJ J £ was 
taken to be 10. We used the damping coefficient a = 0.04 
and the ac bias frequency lu — 0.03 as these values are 
very close to the experimental ones (see the next sec- 
tion). The obtained ip{x,t) profile was used to find the 
voltage (proportional to the fluxon velocity u) across the 
junction. The averaging of velocity was performed over 
the time interval which is a multiple of the period of the 
ac frequency. We used two programs with different nu- 
merical schemes and averaging algorithms, and obtained 
identical results. 




fluxon velocity u 

FIG. 2. Numerically simulated dc driving force vs velocity 
7 (it) curves for various amplitudes of ac drive e indicated in 
the legend. For comparison, dashed line shows the theoretical 
curve for e = 0.1 taken from Fig. Q 

The 7-u curves for different values of the ac bias am- 
plitude £ are shown in Fig. ||. The main effect of the ac 
induced damping which depends on the fluxon velocity 
is clearly present in numerical data. 

Moreover, as the ac bias is small (e < 0.2) the qualita- 
tive agreement between analytical results obtained above 
using the perturbation theory and direct numerical sim- 
ulation is rather good practically in a whole region of 
fluxon velocities u < 1 (see Fig. |^). However, as e in- 
creases the theoretical analysis overestimate the ac in- 
duced damping in the region of small u and underestimate 
it in the opposite regime of large u. Most probably, the 
reason for that is that we neglect the interaction of the 
fluxon with the plasma oscillations. The direct numerical 



simulations also show that for high ac bias (e = 0.5, 0.6, 
and 0.7) the state of moving fluxon becomes unstable, 
and the simulated curves in their upper range end with 
a point of switching to high voltages. 

IV. EXPERIMENT 

To observe the predicted effect we carried out real mea- 
surements using Nb/Al-AlOaj/Nb long annular Joseph- 
son junctior£3 with a trapped magnetic fluxon. The 
junctions were underdamped with the typical value of 
a ~ 0.04 at r = 4.2 K. The junction has the mean 
radius i? = 46 /im and width w — 5 /im. The critical cur- 
rent density jc was about 100 A/cm^, which corresponds 
to the normalized junction length £ of about 10. A single 
magnetic fluxon (Josephson vortex) was trapped in the 
junction by applying a small dc bias current during the 
transition from the normal to tjhe_auperconducting state 
while cooling down the sampleB't2l 
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FIG. 3. The experimentally measured I-V curves of an 
annular LJJ with trapped fluxon for different values of mi- 
crowave ac drive of the frequency 1.4 GHz at power levels 
varying from 4.3 dBm to 20.6 dBm. The lowest curve corre- 
sponds to the autonomous I-V curve (zero ac power). 

For the external ac bias the microwave radiation with 
different frequencies in the range from 1 GHz to 20 GHz 
was used. Here we present the results only for / = 
u!/{2tt) — 1.4 GHz which corresponds to the lu ~ 0.03a;p. 
The microwave signal was applied by an antenna placed 
close to the sample, so that ac bias current was induced 
by the antenna in the bias leads of the junction. The 
typical power W of the applied microwave radiation was 
ranging from few pW to several tens of mW, referenced 
to the top flange of the sample holder. 

We measured I-V curves of the annular Josephson 
junction with one trapped magnetic fluxon at different 
ac power levels W. These I-V are shown in Fig. |^. The 
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measured I-V curves qualitatively look very similar to 
the 7(m) dependencies obtained for the fluxon in Sec- 
tions II and III. Clearly, ac induced voltage-dependent 
damping is observed (cf. Figs. |l| and|2|). 
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FIG. 4. The numerically simulated 7-u curves in the pres- 
ence of magnetic field h — 0.05. The amplitude of ac drive 
was chosen e = 0.3 and 0.4, and other parameters are the 
same as in Fig. hi Arrows identify Shapiro steps. 



VI. CONCLUSIONS 

We have reported a theoretical and experimental study 
of the fluxon motion in the presence of dc and ac bias. By 
making use of the analysis and the direct numerical sim- 
ulations we obtain the ac induced nonlinear damping of 
a fluxon. We also measured the current-voltage charac- 
teristics of an annular LJJ with the trapped fluxon. We 
observed a such ac induced damping as a large ascent of 
the I-V curve in the presence of externally applied mi- 
crowave radiation. Because the response to the ac bias 
increases as the intrinsic damping parameter a reduces 
(Eq. (H) and Fig. 2), the observed effect can be use- 
ful as a method to detect weak microwave radiation. In 
addition, this effects allows to accurately measure the 
damping parameter a when it gets extremely small at 
low temperatures. 
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V. DISCUSSION 

The /- V curves shown in Fig. ^ display many small 
steps which we explain by the presence of small non- 
uniformity of the LJJ. In fact these steps are nothing 
else but conventional Shapiro steps induced by ac- drive. 
The imperfection of the junction can be due to either 
self-field effects, infiuence of the dc bias leads, or small 
technological inhomogeneities. 

To check the influence of non-uniformity on simulated 
characteristics, we performed additional numerical simu- 
lations with a tiny external magnetic field h applied in the 
junction plane. It is a well known that a such magnetic 
field GCjeates the coordinate dependent potential for the 
fluxonETEl, and therefore breaks the translational symme- 
try of the system. This leads to the change of 7-u curve 
such that the fluxon rotation frequency in the annular 
system gets locked to the frequency of ac drive in some 
ranges of 7. Physically, this effect corresponds to the 
appearance of ac-induced voltage locking steps known as 
Shapiro steps for Josephson junctions. 

Figure || clearly shows that in the presence of magnetic 
field there are numerous Shapiro steps of a small mag- 
nitude that appear in the simulated 7-w curve. Yet the 
pronounced effect of ac induced fluxon damping remains 
very evident. 
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